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Abstract. - We develop a Chern-Simons theory to describe a two-dimensional electron gas in 
intermediate magnetic fields. Within this approach, inhomogeneous states emerge in analogy 
to the intermediate state of a superconductor. At half filling of the highest Landau level we 
find unidirectional charge- density- wave (CDW) solutions. With a semiclassical calculation we 
give an intuitive explanation of the change of CDW orientation in the presence of an in-plane 
magnetic field. An anisotropy in the electron band mass is suggested as a possible source of the 
reproducible orientation of the CDW. 



In two-dimensional electron systems both interactions and disorder have surprising and 
striking consequences. A strong perpendicular magnetic field quenches the kinetic energy and 
allows for the observation of both fractional quantum Hall effect and composite fermions in 
the lowest Landau level (LLL), where correlations among electrons are especially important. 
In higher Landau levels on the other hand, interest has mainly been devoted to reentrant 
localization-delocalization transitions causing the integer quantum Hall effect, with interac- 
tions thought of as modifying its non-universal features. The discovery of a pronounced 
resistivity anisotropy close to half filling of higher LLs [El || with filling factor v > 4 has 
recently attracted attention to interaction effects in this regime. The transport anisotropy is 
consistent with the formation of a unidirectional CDW state which was predicted theoretically 
on the basis of Hartree-Fock calculations || ^, || . Numerical exact diagonalization studies also 
support the idea of a CDW formation . The orientation of this CDW structure is related to 
the anisotropy observed in resistivity [Q* || . 

In this letter, we show that a CDW ground state arises naturally in the framework of a 
Chern-Simons (CS) theory The CS approach successfully describes the fractional quantum 
Hall effect and composite fermions in the LLL. In intermediate LLs it is complementary to 
Hartree-Fock calculations which become exact in the limit of high Landau levels ||. Our 
bosonic CS theory has solutions similar to the intermediate state of superconductors and 
allows for an intuitive interpretation of the striped phase: in the highest, partially filled LL 
(HLL) domain walls separate incompressible "superconducting" stripes with local filling factor 
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= 1 from empty stripes (^hll = 0) penetrated by flux. The CDW wave vector is 
determined by the competition of the domain-wall energy favoring a long wave length and the 
Hartree energy which generally increases with CDW wave length but has minima for special 
wave vectors due to the modulation of the electron charge density by the wave-function form 
factor. The CDW wave length and structure are calculated numerically within a mean-field 
approximation. Furthermore, we calculate the wave-function form factor in the presence of an 
in-plane magnetic field using a semiclassical approximation. We find that for a thin quantum 
well the in-plane component of the field is parallel to the axis of high resistivity, in agreement 
with experiments |ll| and numerical calculations fl2"| , 13 . The band mass anisotropy is 



shown to provide a mechanism determining the the experimentally reproducible orientation of 
the CDW wave vector for perpendicular fields. 

Chern-Simons theory. — To employ the CS approach we first map the problem of a partially 
filled iVth LL onto an effective LLL problem. The scattering of electrons between different 
LLs can be neglected when the Bohr radius ob = (47rfi, 2 eoe r )/(me 2 ) is much smaller than the 
magnetic length I — ^Jh/eB, a condition satisfied in experimentally investigated systems. We 
describe screening by the completely filled LLs by a wave- vector dependent dielectric constant 
— e o£ r {l + (2/gaB)[l — J§(<7i? c )]} @ with J denoting the zeroth order Bessel function 
and R c = \/2N + 1 £ the cyclotron radius. The intra-LL matrix elements of the Coulomb 
interaction depend on the LL index N only via the Lagucrre polynomial L^ in the form factor 
of the single particle eigenstates |l5| . One obtains the same matrix elements for states in the 
zeroth instead of the N-th LL but with the effective interaction V cS (q) = ^^L 2 N {Pq 2 /2). 
The Laguerre polynomial can be interpreted as the quotient of the form factor of the iV-th 
and the zeroth LL as Ljv = Fn/Fq. Since all LLs span equivalent Hilbcrt spaces, the new 
interaction V cR in the LLL gives a faithful representation of the original problem. 

The Fermi-Dirac statistics of electrons is mimicked by bosons carrying one flux quantum 
$o = h/e with them. The flux quanta are attached by a singular gauge transformation || with 
a statistical gauge field a(r) tied to the density of bosons p(r) = (jy(r)4>(r) via the condition 
V x a(r) = $op(r). The bosonic CS Hamiltonian is given by 

H = { d 2 r (j>Ur)— (-ihV + e5af <f>(r) + I [ d 2 rd 2 r' 5p(r)V cS (r - r')5p(r') , (1) 
J 2m 2 J 

with <5a = A + a, 5p = <3>^<I> — Itjll-B/^o, and e the charge of the proton. 

Mean-field analysis. — We perform a mean-field (MF) analysis of the bosonic CS problem 
and replace the operators by a complex field. The kinetic energy (of the bosons, not to be 
confused with the quenched kinetic energy of the electrons) alone would be minimized by 
the homogeneous saddle points <\> = y/B/$>o, 6& = and <fr = 0, da = A corresponding to a 
completely filled or empty LL, respectively. Without the interaction potential the ground state 
of a half filled LL would clearly be one domain with filling factor ^hll = 1 and another domain 
with ^hll = 0. However, due to the long-range nature of the Coulomb interaction, these large 
domains break up into smaller ones in such a way that the sum of domain-wall and Coulomb 
energy is minimal, and a charge-density wave is formed in analogy to the intermediate state of 
a superconductor. The effective field V x Sa is expelled from the "superconducting" regions 
with f hll = 1 ■ For average fillings I^hll = 1 — e close to one on the other hand, the solution 
with lowest energy is given by a flux-tube array in analogy to phases of type-I superconductors 
in the intermediate state, where each flux tube can carry a single or multiple flux quanta. For 
the electron problem this is equivalent to the formation of a Wigner crystal. 

To make the physics of the mean-field solution transparent, we express the wave function by 
amplitude and phase, introduce a gauge invariant "velocity" Q, and express the boson current 
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density j in these new variables 

0(r) = \/J /(rKe(r) ' Q( r )=^(ve + ^<5a^ iW^fQW- (2) 

We enforce the constraint on a(r) with the help of a Lagrange multiplier A (which is the 
rescaled zero component of the CS field) and find the saddle-point equations 

-/" + ^Q 2 / + 2(«h-A)/ = , Vx Q = {v + f-\)i , VA = l/ 2 Qxz (3) 

with the vortex density v(r) = £ 2 V x V0(r). The Hartree potential «h is expressed in 
units of £ 2 frui c . One immediate consequence of the saddle-point equations is the screening 
of an electrostatic potential by the CS superconductor in analogy to the screening of a 
vector potential by a conventional superconductor. A negative potential causes a reduction of 
i^HLL = f 2 from one and induces a current via the second part of Eq. (||). This argument can 
be made quantitative by describing the current response to an electric field with the help of 
the Hall conductivity cth = e 2 /h as j = ohE x z. From the third part of Eq. (|) one finds 
VA = (e/fiw c ^ 2 )E leading to an exact cancellation of the external potential in the Schrodinger 
equation for the superfluid density f 2 . Via this mechanism a Bose condensation in a spatially 
varying external potential is possible. Note that the bosonic wave functions are not restricted 
to the LLL in the variational approach. This is an valid approximation in the limit where 
the electron-electron interaction is weak compared to the inter-Landau-level spacing fli| and 
a mixing of Landau levels is negligible. The same assumption is used also in the Hartree-Fock 
approaches [Q, ^) in order to justify that the ground state does not involve single-particle 
states above the partially filled highest LL. Since on the scale t of the typical particle distance 
the strength of the bosonic interaction V is comparable to that of the original fermionic 
interaction (for N not too large as in experiments), the ground state of the bosonic problem 
(Q) will lie to a good approximation in the LLL, even if we search in an extended space of 
functions including higher LLs. Dropping the variational constraint in the bosonic problem 
therefore induces inaccuracies of the same order of magnitude as ignoring the Landau-level 
mixing in the fermionic problem. As can be seen from the variational result below (Fig. [l]), the 
local filling factor overshoots only by a few percent over one, which indicates that the solution 
lies to a very good approximation within the LLL although the constraint is not imposed. 

ID numerical solution. — In order to confirm the physical picture developed above and 
to provide a comparison of our CS approach with previous Hartree-Fock approaches [Q, [|, 
we have numerically determined the MF solution at half filling (Phll = §)• This was 
done approximately by assuming a unidirectional modulation, i.e., that the charge density is 
modulated with a period A in £ direction and is constant in y direction. In this approximation, 
vortices cannot be taken into account as point-like objects. Instead, all vortices within a period 
have to be concentrated on lines, v(r) = —S(x — (n + 1/2)A). For symmetry reasons, 
currents then flow only in y direction, Q(r) = Q y (x)e y . To obtain the MF solution, we use 
a variational ansatz f{x) — ^ fe cos(27rfca;/A) with the constraint f((n+ ^)A) = at the 
vortex positions. In a first step, we minimize Ti for given A by varying the coefficients fk with 
the constraints of half filling and d x Q y {x) = v + f 2 — 1. Then the resulting energy H(A) is 
minimized as a function of A to obtain the CDW period. The choice of ob = £/ 1 \fv fixes the 
interaction strength and is representative for experiments Jl], |l^] . A full 2D numerical solution 
and the analysis of gauge field fluctuations || is beyond the scope of the present analysis. 

We present our results for filling factors v — 2N + \ in the Landau levels N = 8 and 
N = 2. The structure of the MF solutions is depicted in Fig. |l|. For N = 8, f 2 shows a clear 
separation between full (f 2 w 1) and empty (f 2 sa 0) regions. The width of the transition 
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Fig. 1. - Structure of the MF solution: local filling factor of bosons f 2 (full line), local filling factor 
of electrons p e (dashed line; the charge density is — ep e B/$o), and charge-current density j e (dotted 
line; in arbitrary units) as function of x/l (a) for N = 8 with the optimal period A » 3.1i? c ~ 12. 81 
and (b) for N = 2 with the optimal period A « 3.9R C ps 8.7£. 



between the full and empty region is of order £. Therefore, the separation between full and 
empty regions is blurred out for smaller N (for N = 2, see Fig. |l|b). The actual charge and 
charge-current density, which are obtained after convoluting / 2 and f 2 Q y with the appropriate 
form factors are much smoother than / 2 and vary only on the scale of R c . The charge density 
within this LL shows a relative variation of about 20% for N = 8. This magnitude is roughly 
consistent with Hartree-Fock calculations j|, || for large N. For N — 2 the charge density 
shows a stronger modulation of about 60%. Thus, although charge and charge-current density 
are modulated only on the scale A for all N, the relative variation of charge density is found 
to increase with decreasing N. 

The energy per electron as a function of the periodicity is displayed in Fig. ||. The energy 
scale in our saddle-point approximation is the cyclotron energy. We expect that gauge-field 
fluctuations will renormalize it to the Coulomb energy e 2 /47reoe r €. The optimal periodicity 
Ajv, given by the global minimum of this function, is found at As ~ 3.1i? c and A2 ~ 3.9i? c . 
Thus, the periodicity is larger than expected from the first zero of the form factor at A w 2.6i? c 
in agreement with HF calculations j|, || . However, we find that An / R c increases notably with 
decreasing N. In our approach this effect can be understood from the increasing weight of the 
"kinetic" or domain-wall contribution within the CS energy. 

Tilted magnetic field. — So far, we have studied the formation of a CDW by electrons in 
high Landau levels within the CS approach. In experiments |Q, ^| the longitudinal resistivity 
appears to be systematically higher in [110] directions as compared to [110] directions. This 
raises the question of what mechanism breaks the fourfold crystal symmetry and orienting the 
CDW wave vector parallel to the [110] directions. For this reason resistivity was studied in tilted 
magnetic fields ]lo[ |ll| and a sufficiently strong in-plane component B|| of the magnetic field 
was found to orient the directions of large resistivity parallel to B||. We present a mechanism 
for the coupling between the CDW wave vector and B|| that involves the dynamics of the 
electrons in the z direction perpendicular to the plane. We model the confinement of the 
electrons in this direction by a harmonic potential with eigenfrequency loq. While similar 
models were studied recently in the framework of a Hartree-Fock approach [[l2[ [l3| and of 
Laughlin-like states fl6| , we use a simple semiclassical approach to obtain the anisotropic form 
factor of cyclotron orbits in tilted fields, which then can easily be combined with the theory 
of untiltcd fields. 

The classical equations of motion for a single electron subject to a tilted magnetic field 
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Fig. 2. - (a) Plot of the energy per electron in units of hu c as a function of A/_R C for N = 8 (full line) 
and N = 2 (dashed line), (b) Plot of an orbit for a field tilted by 45° with respect to the xy plane 
(with TV = 2, uj x = ui z and u>o = 3uj z over a time of about 16T_). Lengths are in units of R c . The 
cyclotron mode has its long axis parallel to the in-plane field (in x direction), whereas the oscillator 
mode has its longer axis in y direction. 



B = (_B||,0,.B z ) and a harmonic confining potential are linear and possess two eigenmodes 

with frequencies lo\ = \{oJq + w\ + L)%) ± gVA^o + ^1 + w k) 2 ^o^zi where and w z are 
the cyclotron frequencies associated with the in-plane and out-of-plane components of the 
magnetic field, which typically are small compared to luq. For vanishing Bu, the modes with 
frequency uj+ « u)q[1 + oj^/(wq — or u>it f» — w^/(wq — w^)] become a pure oscillation 
in 2 direction or a cyclotron orbit in the xy plane, respectively. We determine the semiclassical 
orbits from the Bohr quantization condition <f p a ■ r a — 2nh(^ + n a ), where a = ± and we 
choose the ground state, n + = 0, for the oscillator-like mode and the and the Nth Landau 
level, n_ = TV, for the cyclotron- like mode. This choice is unique for ujq > uj z . Subsequently, 
we consider only the projection of the orbits onto the xy plane, which are ellipses with half 
axes of lengths 



n 



R 2 xa = (1 + 2N cr) J— 2 A/'/ 2 2^2 ' R l<? _ ~l R l<T ■ ( 4 ) 

The cyclotron-like ellipse is larger than the oscillator- like ellipse (R a - > R a +)- In addition, 
the cyclotron-like ellipse is longer in x direction than in y direction {R x - > Ry-), whereas 
for the oscillator- like ellipse the opposite holds true (cf. Fig. [|b). Because of this contrary 
deformation the resulting anisotropy of the form factor is delicate. The form factor F(q) is the 
Fourier transform of the density p(r) of finding the electron at a position r = r (i) + r + (t) . It 
factorizes into two contributions from the normal modes, -F(q) = i 7 '_(q)F + (q), which are the 
Fourier transforms of the density p±(r) = ^- J Q ± dtS(r — r±(i)) of the normal mode orbits, 
where T± — 2h/lo± is the orbit period. 

In the presence of the tilted field the form factor is anisotropic and we now discuss the 
location of its first zero at small momenta, which essentially determines the CDW periodicity 
and energy. At small momenta, the zeros of F are those of F_ since the cyclotron-like orbits 
have a larger radius. If one estimates the CDW periodicity from the location of the zeros 
of F- (q) one finds the periodicity proportional to the radius of the cyclotron ellipse in the 
direction of the CDW wave vector. The total energy of the CDW is then minimized if its wave 
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vector points parallel to the long direction of the cyclotron ellipse, i.e., parallel to B||. This 
conclusion is in agreement with experiments |l0| , jOJ . 

Finally, we estimate the energy difference between an orientation of the CDW wave vector 
parallel and perpendicular to the in-plane field in the LL N = 2 (assuming again uj x < u) z <C 
luq). The energy change is related to a change of the CDW period which can be approximated 
by SA/A w (R x -/R c — 1) w (uo x /ojq) 2 . The change of the energy per particle with A is 
estimated numerically as SE « —0.1 hu c SA/A from the domain-wall energy. The tilt energy 
per particle is then given by E t ut = tan 2 (#) (u> z /ljq) 3K ^From the threshold angle 9 w 10° 
[ p] where the CDW orientation follows the in-plane component of the field and a typical ratio 
oj z r*> 0.15 wo we find Etut ~ lmK. Consequently, the intrinsic symmetry breaking mechanism 
in the absence of a tilt must lead to an energy gain of the same magnitude. We now examine 
two possible mechanisms. 

Anisotropy. — In the normal state the electron gas has anisotropic resistivity with "hard" 
axis parallel to the [110] crystallographic direction, i.e., perpendicular to the hard axis observed 
in the CDW state. This normal anisotropy, which has been reported for some of the samples 
studied in M, can be ascribed to scattering processes caused by an interface roughness with 
anisotropic correlation lengths £[iIol > £[no] 10 • The observed surface morphology |L8| is 
consistent with such an anisotropic interface roughness. Since the anisotropy in the CDW state 
is determined by collective pinning, whereas it is determined by single-electron scattering in 
the normal state, the direction of the hard axis is not necessarily identical and therefore 
we now examine the collective pinning mechanism. Variations of the interface height d 
result in an electrostatic potential </>d(q) = \end(q) / e(q) . Here n is the density of positive 
donors and the height fluctuations are assumed to be Gaussian correlated, [|d(q)| 2 ] av = 
ZirdQ^xZy exp[— (q 2 £, x + t)< 2 £ 2 )/2] with correlation lengths £ x , £ y A Q and a typical height 
do ~ 2.8A of the order of a GaAs monolayer. Pinning leads to an energy gain via local 
displacements of the CDW profile which are governed by a smectic elastic response Q] . From 
the curvature of the CDW energy as a function of the wave length we estimate the compression 
modulus K x « 10 11 eVm -2 and the bending modulus K y as (A/4ir) 2 K x . For parameters 
representing the experiments |^] we estimated the pinning energy density in a linear response 
calculation fl9|| , 

/n7re\ 2 f d 2 q [\Mlx + 27r/A,^)| 2 ] av 
^ oughW \2vA) J (2tt) 2 K x q 2 +K y q y ■ (5) 

Note that the pinning strength is given by the disorder correlator near the wave vector 2ir/ A 
and therefore it is exponentially small, [\<j>(q x + 2tt/A, q y )\ 2 ] av oc e -2 ^*' A ) . For £ x 3> A this 
exponential dependence suppresses the pinning energy per particle by orders of magnitude 
below lmK per particle as determined above from the threshold tilt angle. Therefore interface 
roughness is irrelevant for the CDW orientation. In addition it would imply that the CDW 
wave vector is aligned in the direction of the shorter disorder correlation length, i.e., the [110] 
direction ||l7| , |l8| in contradiction to experiments. 

As pointed out in Ref. ]2(J the electric field E = en/2eoe r in [001] direction between the 
electron system and the donor layer provides a different symmetry breaking mechanism. We 
now show that this field generates an anisotropic electronic band mass which induces the 
observed CDW orientation by an energy gain on the appropriate scale. Since Ga and As 
carry opposite partial charges they are displaced such that the GaAs bonds in [11-] direction 
are stretched whereas the bonds in [11-] direction are shortened. From pressure experiments 
[ pT| it is known that a shortening of bonds leads to an increase of the effective band mass 
We estimate the effective mass changes m fi jni /fiioi/ 771 * — •"■ ~ ±10~ 4 by identifying the 



m 



electric field with an effective pressure m2\. This anisotropic band mass affects the CDW 
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via the form factor. A semiclassical analysis implies elliptic cyclotron orbits with half axes 
R[ilo]mia]/Rc — 1 ~ ±10~ 4 which leads to an increased CDW period <JApxo]/[no]/A ~ ±10~ 4 . 
In analogy to the analysis of the in-plane field this implies that the mass anisotropy prefers the 
[110] orientation over the [110] orientation by about lmK per electron in agreement with the 
experimentally observed orientation and the energy scale derived from the tilt experiments. 

Conclusion. - We have developed a bosonic Chern-Simons theory to investigate the 
formation of a CDW state in intermediate Landau levels and the influence of an anisotropic 
interface potential on the CDW orientation. In a semiclassical calculation we have found that 
the CDW wave vector is parallel to an in-plane magnetic field for quantum wells with a hard 
confining potential. The band mass anisotropy is a possible origin for the CDW orientation in 
perpendicular fields. 
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